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Abstract. We report the angular diameter measurement obtained with the VINCI/VLTI instrument on the nearby 
^ I star Procyon A (a CMi A, F5IV-V), at a relative precision of ±0.9 %. We obtain a uniform disk angular diameter 

in the K band of 9\ju — 5.376 ± 0.047 mas and a limb darkened value of &ld ~ 5.448 ± 0.053 mas. Together with 
the HIPPARCOS parallax, this gives a linear diameter of 2.048 ± 0.025 Dq. We use this result in combination with 
Q,^ , spectroscopic, photometric and asteroseismic constraints to model this star with the CESAM code. One set of 

• modeling parameters that reproduces the observations within their error bars are an age of 2 314 Myr, an initial 

^■f^ ' helium mass fraction Yi — 0.301 and an initial mass ratio of heavy elements to hydrogen — 0.0314. We also 

. computed the adiabatic oscillation spectrum of our model of Procyon A, giving a mean large frequency separation 

of Afo ~ 54.7 /iHz. This value is in agreement with the seismic observations by Martic et al. (119991 l?001t . The 
interferometric diameter and the asteroseismic large frequency spacing together suggest a mass closer to 1.4 Mq 
than to 1.5 Mq. We conclude that Procyon is currently ending its life on the main sequence, as its luminosity 
O ' class indicates 

w 
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1. Introduction Matching Procyon's position in the HR diagram has 

. : been recognized as of great difficulty by Guenther & 

^ Procyon A (a CMi, HD 61421, HIP 37279) is among the Oemarque CM3 hereafter GD93). The reason is the 

■ ■ ' brightest stars m the sky and is easily visible to the p^^^jy ^^^^^ ^^^^ metallicity Z of the star. We note 

naked eye. This made it an ideal target for a number ^^^^ ^^^^^ computed models only their model b is 

of spectro-photometric cahbration works. It is also a vi- ^j^^^ ^^^^ ^^^^^ frequency spacing measured by 

sual binary star (ADS 6251A) classified F5 IV-V, with M99-M01. Because numerous new studies and observa- 

a white dwarf (WD) companion orbiting the mam com- ^-^^^^ constraints (like the direct diameter) exist today, 

ponent m 40 years. The influence of this massive com- re-examine the status of Procyon A in helium content 

panion on the apparent motion of Procyon was discov- y and in age. We first review and present the adopted fun- 

ered by Bessel iOM^- Guard et al. dUnHD have measured damental parameters of Procyon A in SectlS part of them 

precisely the orbit of the pair, and obtained masses of ^^-^^ ^^^^ j^^, ^1^^ limb-darkening in Sect. where we 

1.497 ± 0.037 Mq and 0.602 ± 0.015 Mq, respectively for ^^^^^^^ interferometric observations and the asso- 

Procyon A and B. It has also beenanasteroseism target ^^^^^^ ^^^^ processing. In Sect.H we detail several models 

since a decade and Martic et al. Ml and 120011 hereafter p^cyon computed with the CESAM code (MorellMJ • 

M99 and MOl) have measured a large frequency spacing of ^j^^^^ ^^^^j^ constrained using the spectroscopic ef- 

respectively 55 and 54 ^Hz. These asteroseismic observa- ^^^^-^^ temperature and the linear diameter value that 

tions provide strong constraints on stellar interior models, ^^^-^^^ ^^^^^ VINCI/VLTI observations. Thus, we 

and on macroscopic stellar parameters like mass and ra- ^^^-^ luminosity for which bolometric corrections 

dius. Comparing the direct interferometric measurements ^^^^^^ uncertain for stellar luminosity class IV-V. A 

of this star to its model diameter is therefore important ^^^^^ asteroseismic frequencies is then presented in 

to cross-vahdate both approaches. gg^^_ ^ ^^^^^^ conclusion. 
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2. Fundamental parameters 

Physical properties of Procyon A are well known thanks to 
the carefully measured orbit of the system and accurate 
HIPPARCOS parallax. Its mass has been recently measured 
by Girard et al. ( WW) at Ma = 1.497 ± 0.037 Mq, adopt- 
ing a parallax of 283.2 ± 1.5 mas. This mass is computed 
using a large number of observations of inhomogeneous 
quality. The same authors when they limit their obser- 
vation sample to the excellent images obtained with the 
WFPC2 found Ma = 1.465 ± 0.041 Mq. If we replace the 
parallax used by Girard et al. H2000|l by the parallax from 
the HIPPARCOS catalogue, 285.93 ± 0.88 mas, the sum of 
the masses of the binary decreases by 2.9 %. Consequently, 
masses become respectively Ma = 1.42 ± 0.04 M© for 
Procyon A and Mb = 0.575 ± 0.017 Mq for the white 
dwarf Procyon B. AUende Prieto et al. H2002I hereafter 
AP02), using the HIPPARCOS parallax of Procyon A, have 
estimated its mass to be 1.42±0.06 Mq. From these works, 
we can conclude that the mass of Procyon is probably be- 
tween 1.4 and 1.5 Mq. Consequently, we investigate these 
two values in our study. The other adopted stellar param- 
eters are summarized in Tabled 

The bolometric luminosity is given by Steffen (19 8 51 
log(L*/L0) = 0.85 ± 0.06. The photospheric properties 
of Procyon have been carefully studied by AP02, who de- 
rived an effective temperature Toff — 6512 K. This value is 
based on the angular diameter measured in the visible by 
the Mark III optical interferometer (Mozurkewich et al. 
11991 1. Detailed 3D stellar atmosphere simulations have 
led AP02 to correct this value of Toff by 80 K. Following 
their results, we adopt Tcff = 6530 ± 50K in our com- 
putations with a surface gravity of \ogg = 3.96 ± 0.02. 
It is interesting to notice that over the last twenty years, 
the effective temperature estimates for Procyon are spread 
between 6545 and 6811 K (Cayrel et al. ll997|l . Adopting a 
color index of B — V = 0.421 and using the calibration by 
Alonso H1996|l we find T^b = 6551 K which is very close to 
the value proposed by AP02. 

The iron abundance at surface of [Fe/H] = — 0.05±0.03 
dex with respect to the solar one is also taken from AP02. 
Therefore we adopt for the modeling of the internal struc- 
ture of Procyon A a chemical mixture which is calculated 
from the iron abundance using the approximation: 



[Fe/H] 



log 



(1) 



log(| 



where {^^)^ is the iron mass fraction within Z. We use 
the solar mixture of Grevesse & Noels (|1993(l then (y )q = 
0.0245. As we shall see in Sect. 01 we adopt at the age of 
Procyon A a surface abundance {§) ^ = 0.0217 ± 0.0017. 
Note that this adopted error bar of AP02 is very small. 
Because the white dwarf (WD) Procyon B has already ex- 
perienced the AGB phase, it is not excluded that material 
coming from this star can have contaminated the narrow 



surface convection layer of Procyon A. However, there is no 
evidence of such an effect in the published table of abun- 
dances of Procyon A fSteffen ri985|l . Enriched elements like 
"s" or other elements coming from Procyon B during its 
post-AGB phase do not seem to have migrated as gas or 
dust in the atmosphere of Procyon A. The mass transfer 
during this short period of the WD progenitor life could 
not exceed the accretion of the wind passing close to the 
star. Therefore, except if the two stars have filled in their 
Roche lobe, the mass of Procyon A has only been changed 
by a negligible ammount. 

The projected rotational velocity of Procyon A has 
been determined by many authors who derived values 
smaller than 5.0kms~^. Considering in particular the 
value proposed by AP02, vs'mi = 3.2 ± 0.5km.s~^, we 
conclude that Procyon A is a slow rotator and we ne- 
glect the rotational velocity in our modeling of its internal 
structure. The semi-major axis of the Procyon orbit be- 
ing a ~ 4'.'5 (Girard et al. l2000|l . the distance between the 
two companions amounts to ~ 16 AU. The tidal interaction 
between the two components is therefore negligible. 

3. Interferometric angular diameter 

3.1. Instrumental setup and data processing 

The European Southern Observatory's Very Large 
Telescope Interferometer (Glindemann et al. I2000|l is op- 
erated on top of the Cerro Paranal, in Northern Chile, 
since March 2001. The observations reported here were 
done with two test siderostats (0.35 m aperture) and the 
VINCI beam combiner (Kervella et al.|l)00, 2003a|). The 
Procyon visibility measurements were all obtained on the 
B3-D3 baseline, 24 m in ground length. We used a regular 
K band filter (A — 2.0 — 2.4 /im) for these observations. 
The effective wavelength changes slightly depending on 
the spectral type of the observed target. For Procyon, we 
estimate Aefj = 2.178 ± 0.003 /xm. This error bar adds 
quadratically a relative systematic uncertainty of 0.15% 
on the final limb darkened angular diameter error bar. 

The raw data processing has been achieved using an 
improved version of the standard data reduction pipeline 
(Kervella et al. '2003^ , whose general principle is similar 
to the original FLUOR algorithm (Coude du Foresto et 
al. I1997|l . The two calibrated output interferograms are 
subtracted to remove residual photometric fluctuations. 
Instead of the classical Fourier analysis, we implemented 
a wavelet based time-frequency analysis (Segransan et al. 
I1999|l . The output of the processing pipeline is a single 
value of the squared coherence factor ij? for each series of 
500 interferograms and the bootstrapped error bar. 



3.2. Visibility calibration 

Sirius was chosen as the calibrator for Procyon obser- 
vations as its limb darkened (LD) angular diameter has 
been measured recently with high precision by Kervella 
et al. I|2003d|l . using VINCI and the 66m EO-Gl base- 
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Table 1. Relevant parameters of a CMi (Procyon) and its 
calibrator aCMa (Sirius). For Sirius, see also Kervella et 
al. Il2008dll 





a CMi 


aCMa 
HD 4SQ1 5 


my 


0.34 


-1.47 


rriK 


-0.65 


-1.31 


Spectral Type 


F5IV-V 


AlVm 


MiMe) 


1.42 ±0.04 


2.12 ± 0.02 


Teff (K) 


6530 ± 50 


9900 ± 200 


logs 


3.96 


4.30 


[Fe/H] 


-0.05±0.03 


0.50±0.20 


V sin i 


3.16±0.50 


16.0±1.0 


n"^ (mas) 


285.93 ± 0.88 


379.22 ± 1.58 


6'uD (mas) 


5.376 ±0.047 


5.94 ±0.02*= 


6'ld (mas) 


5.448 ± 0.053 


6.04 ± 0.02*^ 



" Parallax values from hipparcos (Perryman et al. 119971 . 
^ Sirius 6uu and Sld are taken from Kervella et al. Il2003dll . 



line. Table gives the relevant parameters of Procyon 
and Sirius. The interferometric efficiency (IE) of the in- 
strument was estimated from the values obtained on 
Sirius a short time before and after Procyon (typically 15 
minutes) . Thanks to the very high SNR of the Sirius inter- 
ferograms, we were able to obtain very small relative errors 
on the IE values, typically 0.2% for a series of 500 interfer- 
ograms. The error bar on the angular diameter of Sirius 
is treated as a systematic error in the calibration process, 
and is therefore not averaged for multiple observations. 
For the visibility fit, we took into account simultaneously 
the limb darkening and the bandwidth smearing, as de- 
scribed in Kervella et al. I|2003bll2003djl . 

Our February 2003 observation campaign of Procyon 
resulted in a total of 53 calibrated V'^ measurements, from 
23 256 processed interferograms. In the calibration pro- 
cess, we separated clearly the statistical and systematic 
error contributions. This is essential to avoid underesti- 
mating the final error bars of the fit. The latter cor- 
responds to the uncertainty on the angular diameter of 
Sirius. In spite of the fact that the angular size of Sirius 
is larger than Procyon's, the small error on its measured 
value (±0.3%) translates into a small calibration error. 
As a consequence, the global error on the angular size of 
Procyon is dominated by the statistical contribution by a 
factor of two compared to the systematic calibration error: 
the uncertainty on the angular size of Sirius is not limiting 
our final precision. The visibility points as a function of 
the projected baseline are presented on Fig. ^ and their 
residual scattering around the best fit model on Fig. |21 

3.3. Limb darkened angular diameter and linear size 

The intensity profile /(/x)//(l) that we chose for Procyon 
was computed by Claret et al. 1)2000(1 . based on the 
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Fig. 1. Squared visibility measurements obtained on 
Procyon and best fit model (solid line). 
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Fig. 2. Scattering of the Procyon measured visibilities 
around the best fit model. The plotted error bars are 
the quadratic sum of the statistical and systematic con- 
tributions. 

ATLAS9 model atmospheres fKurucz [l992|l . It is a four 
parameters law: 

4 

/(:r)//(l)^l-^a,(l-xt) (2) 

k=l 

which coefficients are ai — 0.5089, a2 = 0.0822, 03 = 
—0.3978, and 04 = 0.1967. The parameter x = cos{9) 
is the cosine of the azimuth of a surface element of the 
star. The resulting profile is shown on Fig.|3| AP02 have 
identified an influence of the convective granulation of 
Procyon on its limb darkening at visible wavelengths. 
Though, the amplitude of this effect is already very small 
at A = 1.0 /im, and we neglect it for our observations at 
A = 2.2 /im. 

We obtain directly the LD angular diameter of 
Procyon from a classical minimization, 6'ld = 5.448 ± 
0.053 mas. Using a simple uniform disk model, we find a 
value of 0UD = 5.376 ± 0.047 mas. 

We can compare the LD value with the previously pub- 
lished interferometric measurements, listed in Table[21 Our 
value of 0LD is compatible with previous measurements, 
that were all obtained using visible wavelength observa- 
tions. It is also in excellent agreement with the average 
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Table 2. Intcrfcrometric measurements of the angular di- 
ameter of Procyon A from the Hterature. 



Instrument 


A (^m) 


6'uD (mas) 


^LD (mas) 


Intensity interf.^ 


0.45 


5.10 ±0.16 


5.41 ±0.17 


Mark III'' 


0.80 


5.26 ±0.05 




Mark 111*= 


0.45 


5.14 ±0.05 




Mark Iir 


0.80 




5.46 ± 0.08 


NPOr 


0.74 


5.19 ±0.04 


5.43 ±0.07 


VINCI/VLTI 


2.18 


5.38 ± 0.05 


5.45 ±0.05 


Average value 






5.445 ± 0.035 



Hanbury Brown et al. II1974II . 
Mozurkewich et al. J1991II . 
" Nordgren et al. ll^nnTt . 



S 80% 



1.0 0.8 0.6 0.4 0.2 0.0 

X 

Fig. 3. Intensity profile of Procyon, from the four param- 
eters law of Claret IpW))! . 

of all published values. From the VINCI/VLTI value of 
^LD and the hipparcos parallax, we deduce a linear di- 
ameter of 2.048 ± 0.025 Dq. Computing the average of 
all published 0ld values, we obtain a linear radius of 
= 2.047 ± 0.020 R© for Procyon A, statistically iden- 
tical to our result. 

4. Modeling 

The modeling of Procyon A is based on the matching 
of spectroscopic effective temperature, surface metallicity 
and intcrfcrometric radius by a "satisfactory" evolution- 
ary model for a given mass. The constraint between i?+ 
and Tcff within their error bars is illustrated in Fig. 0] by 
the hatched parallelogram, while the constraint between 
the Teff and is illustrated by the dashed rectangle. That 
emphazises the advantage of the use of the measured i?* 
instead of the L^, which depends of photometric calibra- 
tions and bolometric corrections. 

Table |2| shows the characteristic and initial modeling 
parameters of the star. The ordinary assumptions of stellar 
modeling are made, i.e. spherical symmetry, no rotation, 
no magnetic field and no mass loss. The relevant nuclear 
reaction rates are taken from the NACRE compilation 
by Angulo et al. H1999fl . The equation of state adopted 
is taken from Eggleton et al. H1973|l , and the opacities are 
from the OPAL database (Iglesias & Rogers [TI)96(l . using 
the Grevesse & Noels p993f) mixture. The microscopic dif- 
fusion is described using the formalism of Burgers H1969II 



with the resistance coefficients of Paquette et al. p98tifl . 
We take into account the radiative diffusivity as recom- 
mended by Morel & Thevenin 1)20021 , that limits the cfli- 
ciency of the microscopic diffusion in outer-layers of stars 
with intermediate masses. The atmosphere is restored us- 
ing Hopf's law (Mihalas ll978|l . The definition of the radius 
of the star is the bolometric radius, where T(i?i) = Tgg- 
In the convection zones the temperature gradient is com- 
puted according to the MLT^m convection theory with a 
mixing length parameter of A = 1 (Canuto & Mazzitelli 
[T99L 1992^. 

Following the discussion of Sect. [3 we investigate the 
sensitivity of our models to a variation of the mass of 
Procyon in the range 1.4 to 1.5 Mq. We obtained a "sat- 
isfactory" model a (Fig. 2)), detailed in Table |31 which 
reaches the hatched parallelogram corresponding to an age 
of 2 314 Myr. Its evolutionary state corresponds closely to 
the end of the main sequence, when the convective core 
has disappeared owing to the exhaustion of the hydrogen 
at center. We present on Fig. [S] the variation for model 
a oi Xs, Ys and of [Fe/H]s as a function of the age. From 
the ZAMS until about 1 Gyr, due to the gravitational 
settling, the surface abundances of helium and heavy el- 
ements (resp. hydrogen) decrease {resp. increases). After 
1 Gyr, the density in the envelope decreases due to the 
enlargement of radius. This leads to an increase of the ra- 
diative mixing, with the consequence of a dredge-up which 
increases {resp. decreases) and [Fe/H]s {resp. Xs) at the 
surface. 

We also computed a "satisfactory" model b without 
microscopic diffusion in order to estimate its effect on the 
age of the star and we found that suppressing the diffusion 
increases the age by ~ 400 Myr. The evolutionary state of 
this model is the beginning of H burning in shell. 

Following the prescriptions of Schaller et al. 1)1992(1 we 
have computed additional models that include overshoot- 
ing of the convective core of radius fi!co over the distance 
Ov = A min(i/p, 7?co)- It has been impossible to find a 
"satisfactory" model reaching the narrow hatched area of 
the HR diagram with parameters A > 0.03. 

Finally, we computed a model c with Girard et al.'s 
(2000) mass : 1.50 Mq and find an age of 1.3 Gyr without 
overshoot. We also computed a model with an overshoot 
of 0.15 in agreement with Ribas, Jordi & Carme l)20()()|l 
for stars of mass lower than 1.5 M© resulting in a model 
which increases the age of the star to 1.55 Gyr. 

The parameters for all models are listed in Table 13 
Fig. 2] shows the corresponding evolutionary tracks in the 
HR diagram. We note that our models fit the value of 
log 5 = 3.96 proposed by AP02. 

5. Asteroseismology 

Asteroseismic observations of Procyon A have been re- 
ported by M99 and MOl and indicate an excess of power 
around 1 milliHertz with a probable frequency spacing of 
54 to 55 fiUz. 
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Fig. 4. Evolutionary tracks in the HR diagram of mod- 
els of Procyon A (see Tabled. Dashed rectangle delimits 
the uncertainty domain for luminosity and effective tem- 
perature, while the hatched area delimits the uncertainty 
domain for effective temperature and the interferometric 
radius. 




500 1000 1500 

Age (Myr) 



Fig. 5. Changes with respect to time for Xg, Yg and 
[Fe/H]s for model a of Procyon. 

The narrow convective envelope of the star we con- 
sider may stochastically excite solar-like oscillations. As 
in the Sun, these oscillations will have quasi equidistant 
frequencies at given degree. 

For all the models, we have computed a set of adiabatic 
frequencies of the normal modes for degrees £ = 0,1, 2, 3. 
The set of frequencies of model a are given in the frequency 
range 450 to 1 600 /xHz in Table H 

The large frequency spacing is defined as the differ- 
ence between frequencies of modes with consecutive radial 
order n: Ai>£{n) = Vn^e — t^n-i,i- In the high frequency 
range, i.e. large radial orders, Avi is almost constant with 
a mean value Avo, strongly related to the mean density 
of the star. The computed frequencies are fitted to the 
following asymptotic relation (Berthomieu et al. 11993(1 : 
i^n,e — >^o + Avi{n — uq) + ae{n — no)^. With no = 21, 
£ = and radial order n between 17 and 26 we obtain an 
estimate of the mean large difference Avq. 



Table 3. Procyon A models (without overshoot, see text) 
lying within uncertainty box in the hr diagram. The sub- 
scripts "i" and "s" respectively refer to initial values and 
surface quantity at present day. "c" refers to the central 
value. 



iVioaei 


a 


D 


c 


m/Me 


1.42 


1.42 


1.50 


Yi 


0.3012 


0.2580 


0.345 


Ys 


0.2209 


0.2580 


0.202 


(#). 


0.03140 


0.0218 


0.0450 


(#). 


0.02157 


0.0218 


0.0220 


diffusion 


yes 


no 


yes 




0.00051 


0.00000 


0.2180 


age (Myr) 


2 314 


2 710 


1300 


Teff (K) 


6524 


6547 


6553 




3.960 


3.967 


3.994 


[Fe/H]i 


-1-0.107 


-0.051 


4-0.264 


[Fe/H]. 


-0.055 


-0.051 


-0.043 


log(L/Lo) 


0.8409 


0.8405 


0.8390 


R/Rq 


2.0649 


2.0495 


2.0420 


Avo (/iHz) 


54.7 


55.4 


56.4 



The mean large spacing of model a is 54.7 /iHz 
which fits well observational predictions of M99 & MOl 
(AfQ ~ 54 — 55 /iHz). For the model c, the mean large 
spacing is found to be 56.4 /tHz, significantly different 
from the observed value. A mass as large as 1.5 Mq 
for Procyon A seems therefore improbable. Adopting 
the VINCI/VLTI radius, and looking at the scaling of 
Kjeldsen & Bedding ipH^ . 

~ 134.9^(£^ i^iHz), (3) 

the most straightforward way to decrease Avo is to de- 
crease the mass of Procyon A. 

Because the luminosity L is in relation with Toff 
and with the radius R through the formula L/Lq — 
{R/RQf{T^s/T%Y we find from Eq.Ethc expression of 
the mean large spacing as a function of L and Toff, 

From this formula it is easy, if needed, to estimate the 
derivatives of the mean large spacing with respect to the 
mass, the effective temperature and the luminosity. To 
estimate the dependance between Aj/qsc and Zi we have 
computed a modified model a with AZi = 0.001 and found 
a variation in the parallelogram error box of : ATefj — 
22K, AL/L = 0.008, and Agg e = 33 Myr. There were no 
significant change with Zi of A 

^osc 03(11 be expected 

from Eq. |21 

6. Discussion & Conclusions 

We have reported in this paper our modelisation of 
Procyon A based on observations by long-baseline inter- 
ferometry and asteroseismology. The use of the measured 
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Table 4. Astcroscismic frequencies (/iHz) of Procyon A 
for model a. The correspond to "g" modes which are 
presents in the frequency domain due to the evolutionary 
state of the star. 



£^0 


£=1 


£^2 


£^3 


n 


485.35 


513.36 


532.31 


551.46 


8 


- 


- 


- 


597.20 


★ 


537.93 


563.43 


585.90 


605.61 


9 


591.09 


615.78 


638.60 


657.72 


10 


643.23 


666.79 


689.80 


708.98 


11 


694.17 


717.44 


741.00 


760.87 


12 


745.32 


769.15 


793.75 


814.22 


13 


798.11 


822.24 


847.62 


868.87 


14 


- 




896.38 




★ 


852.14 


876.76 


904.48 


923.72 


15 


906.98 


931.22 


957.38 


978.12 


16 


961.42 


985.26 


1011.19 


1031.92 


17 








1075.73 


★ 


1015.27 


1038.85 


1065.03 


1086.34 


18 


1069.16 


1092.64 


1119.00 


1140.31 


19 


1123.16 


1146.76 


1173.64 


1195.27 


20 


1177.86 


1201.57 


1228.79 


1250.68 


21 


1233.05 


1256.82 


1284.20 


1306.18 


22 


1288.46 


1312.09 


1339.65 


1361.81 


23 


1343.91 


1367.55 


1395.21 


1417.43 


24 


1399.42 


1422.95 


1450.65 


1472.94 


25 


1454.83 


1478.24 


1506.00 


1528.39 


26 


1510.12 


1533.53 


1561.40 


1583.89 


27 



Further progress on the modeling of Procyon will be 
possible when the accuracy on the flux of the star is im- 
proved to less than 1 %. Waiting for such accuracy, the 
uniqueness of the solutions resulting from computed mod- 
els fitting a narrow box in the HR diagram will come from 
future detailed asteroseismic studies. For example, other 
linear combination of frequencies like the small frequency 
spacings (see e.g. Gough [T991|) will constrain the age and 
the mass of the star. We therefore encourage observers to 
progress on the asteroseismology of this star. 

Large uncertainties also come from the adopted chem- 
ical abundance mixture Zg which is still rather uncertain. 
Only a few chemical element abundances are measured 
today and most of them with a low accuracy. This uncer- 
tainty on Zs is the largest source of error on the estimated 
initial helium content Yi and on the age of Procyon. Thus, 
we recommend that surface abundances should be derived 
from 3D atmosphere studies, in particular for oxygen and 
other important donors of electrons. 
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stellar diameter allows to reduce significantly the error bar 
on the luminosity of the star, as it does not require any 
bolometric corrections. This advantage has been demon- 
strated in this study and in the modeling of Sirius by 
Kervella et al. H2003d|l . 

Our model shows that using the given set of param- 
eters and given physics, Procyon A is currently finishing 
to burn its central hydrogen, and is at the phase where 
the convective core is disappearing. The error on the mea- 
sured radius gives a narrow uncertainty of 10 Myr on the 
deduced age. Provencal et al. H2()()2|l have discussed the 
cooling time of the WD Procyon B and found that the 
progenitor ended its lifetime 1.7 ± 0.1 Gyr ago. We derive 
an age of 2 314 Myr for Procyon A. Subtracting the cool- 
ing age of the WD companion to our determination of the 
age of Procyon A leads to a lifetime of 614 Myr for the 
progenitor of Procyon B. This indicates that the mass of 
the progenitor is approximately 2.5 Mq. This value yields 
a mass of ~ O.57M0 for the core of the corresponding 
Thermal-Pulsating-AGB star (for Z=0.020) (Bressan et 
al. I1993|l . which is the minimum possible value for the fi- 
nal mass of the WD (see also Jeffries Il997|l . This estimate 
of 0.57 M0 agrees very well with the mass of Procyon B 
that we deduced from Girard et al. H2000|l . We note that 
the age obtained with model c is shorter than the cooling 
age of Procyon B. This argument suggests a mass lower 
than 1.5 Mq for Procyon A and strengthens the astero- 
seismology results. 
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